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PULSE  COMPRESSION  DEGRADATION  DUE  TO  OPEN 
LOOP  ADAPTIVE  CANCELLATION,  PART  I 


1.  INTRODUCTION 

An  adaptive  canceller  [1]  adjusts  the  auxiliary  channel  settings  so  as  to  null  out  interference  in 
the  main  channel  (Fig.  1).  The  weights  of  the  auxiliary  channels  are  adjusted  to  minimize  the  output 
noise  power  residue.  If  the  noise  environment  is  not  known  a  priori,  these  weight  settings  must  be 
estimated  from  a  finite  set  of  incoming  data  on  the  input  channels.  Thus,  the  weight  settings  will 
have  perturbations  about  the  quiescent  optimum  weight  settings.  These  perturbations  cause,  for  exam¬ 
ple,  the  output  noise  power  residue  to  rise  [1-4],  and  if  the  adaptive  canceller  is  used  in  the  sidelobe 
antenna  canceller  configuration,  the  adaptive  array  antenna  sidelobe  level  to  increase  [5]. 


MAIN  ANTENNA  AUXILIARY  ANTENNAS 


A  pulse  compressor  is  essentially  the  matched  filter  for  a  given  radar  waveform.  Most  of  the 
energy  in  the  received  radar  waveform  is  compressed  into  a  given  single-range  cell  and,  thus,  the  sig¬ 
nal  level  can  be  increased  significantly  for  detection  purposes.  However,  some  energy  does  leak  into 
the  sidelobes  of  the  compressed  pulse  response,  resulting  in  low  gain  in  range  cells  outside  of  the 
given  range  cell.  If  a  target  or  piece  of  clutter  is  large  enough,  it  can  break  through  and  be  detected 
in  these  range  sidelobes,  falsely  indicating  a  target  detection  or  masking  a  real  target.  Thus,  it  is 
highly  desirable  to  maintain  a  low  sidelobe  respon.se. 
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For  the  radar  designer  there  is  the  choice  of  where  to  put  the  pulse  compressor:  before  or  after 
the  canceller.  A  disadvantage  of  placing  it  before  the  canceller  is  that  a  pulse  compressor  must  be 
placed  in  each  antenna  channel  (main  and  auxiliaries)  to  maintain  channel  match.  Another  disadvan¬ 
tage  is  that  the  pulse  compressor  must  have  the  dynamic  range  of  the  interference  (possibly  clutter  and 
jamming)  which  has  yet  to  be  cancelled.  These  disadvantages  do  not  exist  if  the  pulse  compressor  is 
placed  after  the  canceller.  However,  a  disadvantage  of  placing  the  pulse  compressor  after  the  can¬ 
celler  is  that  the  range  sidelobe  levels  of  the  compressed  pulse  increase  because  a  finite  number  of 
samples  is  used  to  compute  the  canceller  weights.  (We  assume  that  the  canceller  weights  are 
estimated  from  a  block  of  input  data  and  applied  back  onto  the  same  input  data  set.)  Hence,  the 
designed  radar  waveform  pulse  compressor  responses  may  not  be  achieved  because  of  the  interaction 
of  the  canceller  with  the  pulse  compressor.  Because  of  the  higher  sidelobes,  unwanted  targets  may 
break  through  and  be  detected.  Significantly,  as  it  will  be  shown,  the  maximum  sidelobe  levels  that 
result  from  interacting  with  the  canceller  are  independent  of  whether  the  input  waveform  was 
compressed  before  or  after  the  canceller  if  the  uncompressed  pulse  is  completely  contained  within  the 
samples  that  are  used  to  calculate  the  canceller  weights. 

This  report  considers  the  degradations  that  result  .vhen  trying  to  match-filter  or  pulse-compress  a 
desired  radar  waveform  after  it  has  been  processed  through  a  particular  adaptive  canceller  called  the 
Sampled  Matrix  Inversion  (SMI)  algorithm  [2].  The  SMI  algorithm  is  an  open-loop,  rapidly  converg¬ 
ing  adaptive  canceller  implementation  whose  noise  power  convergence  rate  is  independent  of  the 
external  noise  environment.  For  many  years,  it  has  been  considered  a  baseline  for  fast  converging 
adaptive  canceller  algorithms.  The  Gram-Schmidt  (GS)  canceller  [4,  6-8]  is  a  numerically  equivalent 
implementation  (assuming  infinite  accuracy)  of  the  SMI  algorithm  with  excellent  performance  simul¬ 
taneously  in  arithmetic  efficiency,  stability,  and  convergence.  In  addition,  it  is  a  good  analytical  tool 
[4,5]  with  which  to  investigate  the  convergence  properties  of  the  SMI  canceller.  Because  of  the 
SMI’s  complexity,  we  assume  that  it  is  implemented  digitally.  Because  the  pulse  compressor  follows 
the  canceller,  we  assume  that  it  is  also  implemented  digitally. 

This  report  is  laid  out  as  follows.  Section  2  briefly  describes  the  GS  canceller  and  Section  3 
presents  the  signal  model  (for  pulse  compression).  Some  past  results  on  GS  cancellers  are  reviewed 
in  Section  4.  Sections  5  through  7  derive  the  analytical  results  and  Section  8  discusses  them. 

For  some  applications,  the  matched  filter  is  replaced  by  a  filtering  scheme  that  reduces  the  range 
sidelobes  at  the  expense  of  signal  gain  at  the  matched  point.  Thus,  a  mismatched  filter  is  used.  It  is 
shown  that  the  results  derived  for  the  matched  filter  are  applicable  to  the  mismatched  filter. 

2.  THE  GS  CANCELLER 

Consider  the  general  A-input  Gram-Schmidt  canceller  structure  as  shown  in  Fig.  2.  Let  .r^(t). 
.ri(/),  ...,  represent  the  complex  data  signals  in  the  0th,  1st,  ....  A  -  1th  channels, 

respectively.  We  call  the  leftmost  input  x^it)  the  main  channel  and  the  remaining  A  -  1  inputs,  the 
auxiliary  channels.  The  main  channel’s  signal  consists  of  a  desired  signal  plus  additive  noise  (i.e., 
internal  plus  external  noise).  Cancellation  of  the  signals  from  external  interfering  sources  relies  on 
the  correlation  of  simultaneously  received  signals  in  the  main  and  auxiliary  channels.  The  internal 
noises  on  each  channel  are  assumed  to  be  uncorrelated  between  channels.  The  GS  canceller  decorre¬ 
lates  the  auxiliary  inputs  one  at  a  time  from  the  other  inputs  by  using  the  basic  2-input  GS  processor 
as  is  shown  in  Fig.  3.  For  example.  Fig.  2  shows  that  Jt/v_i(t)  is  uncorrelated  with  .v)i^’(/). 
jt'i“’(t) . x^Liit)  in  the  first  level  of  decomposition.  Next,  the  output  channel  that  results  from 
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Fig.  2  —  GS  structure 
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Fig.  3  —  Basic  2-input  GS  canceller 


decorrelating  3r/v_i(0  with  a:/v-2(0.  is  decorrelated  from  the  other  outputs  of  the  first-level  GSs.  The 
decomposition  proceeds  until  a  final  output  channel  is  generated.  If  the  decorrelation  weights  in  each 
of  the  2-input  GSs  are  computed  from  an  infinite  number  of  input  samples,  then  the  main  output  chan¬ 
nel  is  totally  decorrelated  with  the  input:  xiit),  XiiO,  . . .  ,  -^0).  Also,  .r},''' ""’(t), 

n  =0,  1,  . . .  ,  N  -  1  are  pairwise  statistically  uncorrelated  or  orthogonal. 

If  an  infinite  number  of  input  samples  is  not  used  then  the  decorrelation  weights  associated  with 
each  2-input  GS  canceller  are  estimated  by  using  finite  averaging.  One  method  of  processing  data 
through  the  GS  canceller  is  called  concurrent  processing  because  the  weights  are  estimated  from  a 
block  of  input  data  and  applied  back  onto  the  same  input  data  set.  (For  a  discussion  of  other  types  of 
GS  processing,  see  Ref.  4.)  This  type  of  processing  is  often  necessary  to  handle  blinking  interference 
sources. 


We  briefly  describe  the  concurrent  GS  canceller.  Let  jcJ,'”’  represent  the  outputs  of  the  2-input 
GSs  on  the  (m  -  1)  level.  Then  outputs  of  the  2-input  GSs  at  the  mth  level  are  given  by 


yim  +  1)  _  (m) 


W 


im)  Un) 
n 


n  =0,  N  —  m 

m  =  \,  2,  ,  N  -  \ 


1 


(1) 


Note  that  jCn’  =  x„.  The  weight  w*,'"’,  seen  in  Eq.  (1),  is  computed  so  as  to  uncorrelate 
n  =0,  1,...  ,  N  -  m  -  1  with  x^l„.  For  K  input  samples  per  channel,  this  weight  is  estimated 


as 


£  x^j(;’i^*{k)xy{k) 

k=\ _ 

k =1 


(2) 


where  *  denotes  the  complex  conjugate  and  |  •  1  is  the  magnitude.  Here  k  indexes  the  sampled  data. 
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Let  xq  represent  the  additive  noise  in  the  main  channel.  For  this  development,  we  make  the  fol¬ 
lowing  five  assumptions; 

1.  The  Tq,  -Vi,  are  identically  distributed  Gaussian  complex  random  variables 

(r.v.). 

2.  These  same  r.v.s  are  samples  from  stationary  processes  with  zero  mean. 

3.  .r„ I  (/ti)  is  statistically  independent  of  jr„, (it 2)  for  ij  ^  k2- 

4.  The  desired  signals  are  not  present  in  the  auxiliary  channels. 

5.  K  >  N. 

It  is  necessary  to  make  the  last  assumption  because  for  A"  <  N  -  \  the  GS  canceller  is  numeri¬ 
cally  unstable  because  of  the  singularity  of  the  estimated  input  covariance  matrix,  and  for  A'  =  A'  — I 
the  output  of  a  concurrent  GS  canceller  is  always  zero  (4].  In  the  following  discussion,  a  normalized 
A-length  multivariate  complex  circular  Gau.ssian  vector  is  defined  to  have  L  elements,  each  of  which 
has  real  and  imaginary  parts  that  are  independent  Gaussian  r.v.s  with  0  mean  and  variance  equal  to 
1/2  (the  magnitude  variance  is  one).  In  addition,  the  L  elements  are  independent  of  one  another. 

Figure  4  presents  simplified  /V-input  GS  canceller  structures  for  concurrent  processing.  The 
notation  GS/(j^  indicates  that  an  A-input  GS  structure  uses  K  samples  from  each  channel  to  compute 
the  weights  interior  to  the  GS  structure.  The  0th  channel  (the  far  left  channel  in  Fig.  2  or  4)  is 
always  designated  as  the  main  channel  and  the  others  are  called  the  auxiliary  channels  (or  just 
AUXs).  Figure  5  represents  the  GS  structure  with  N  orthogonal  outputs  displayed. 


Fig.  4  —  Representation  of  GS  caneeller  with  N  channels  and 
K  samples  per  channel 
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AUX  CHANNELS 


Fig.  5  —  GS  representation  with  N  output  channels 


3.  SIGNAL  MODEL 

A  network  whose  frequency  response  function  maximizes  the  output  peak  signal-to-mean  noise 
power  (o/iV)  .-i'c  is  called  z  '^a’ched  filt"''  Almost  all  radar  receivers  are  designed  by  using  the 
matched  filter  criteria.  If  h{t)  is  the  impulse  response  function  of  the  matched  radar  receiver.  s{t)  is 
the  transmitted  radar  waveform,  and  the  noise  interference  is  white  and  additive,  then  it  can  be  shown 
[9-11]  that 


h(t)  =  — /) , 


(3) 


where  *  denotes  the  complex  conjugate  operation. 

A  sampled  matched  receiver  design  is  based  on  the  same  principle  of  maximizing  S  /N.  We 
sample  the  transmitted  radar  waveform  at  equal  time  intervals  t.  Let  5 1 ,  •  •  ■  •  the  values  of 

the  sampled  transmitted  waveform  or  code  where  L  is  the  number  of  sampled  points  (Fig.  6).  Set 

S  =  (5l,  52 . 5^)^.  (4) 

where  T  denotes  the  vector  transpose  operation.  The  sampled  receiver  applies  a  weighting  vector 
g  =  (;?i-  1?2, - such  that 


y  =  g^s.  (5) 

It  can  be  shown  that  (S  IN)  is  ma.ximized  when 

g  =  s*.  (6) 
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Fig  6  —  Sample  radar  wavelorm 


U(n) 

n  =  0,1,2.... 


MATCHED 

FILTER 
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V(n) 


Fig.  7  —  Sampled  matched  filler 


Figure  7  shows  the  matched  receiver  structure.  Here  U(n)  is  the  received  time  sampled  input 
sequence  consisting  of  signal  plus  noise.  The  received  input  sequence  is  convolved  with  conjugated 
time  reversed  vector  of  s  which  result.-^  in  an  output  .sequence  V{n).  Mathematically,  this  is  written 

i  - 1 

VOi)  =  +  !).  (7) 

/  -0 

The  sampled  matched  filter  response  of  the  desired  signal  vector  s  of  length  L  is  given  by 


/  -m 

E  t 


/  =  I 


m ' 


r(m)  =  ^ 


E 


I  ^  I 


0  <  m  -  i. 


-  (T  -  I)  <  w  <  0. 


(8) 


This  sequence  is  2Z.  -  I  points  long  and  is  often  called  the  autocorrelation  function  (ACF)  of  .v , . 

S  2 . ^  * 
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If  we  define 


r  =  (r(-(L  -  D),  .  .  .  ,  r(-l).  r(0),  r(l) . riL  -  1))^ 

and 


(9) 


sl  0  0 

Sl  - 1  5/.  0 
^L-1  ^L-\ 


0 

0 

0 


5^  = 

S2 

•  •  •  5/. 

0 

52 

•  •  •  -S/.  -  1 

0 

0 

•G 

••• 

(10) 


0  0  0  •  •  •  5 , 


where  r  is  a  2L  -  1  length  vector,  and  5  is  a  Z-  x  (2L  -  1)  matrix  called  the  ACF  matrix  of  s.  we 
can  show  that 


r  =  S's 


(11) 


where  /  is  the  conjugate  transpose  operation. 

We  assume  for  this  analysis  that  the  GS  canceller  processes  data  in  blwks  of  K  data  samples  per 
channel.  Thus,  the  desired  signal  vector  may  be  spread  across  a  number  of  sample  blocks.  How¬ 
ever.  for  this  analysis,  we  assume  the  signal  vector  of  length  L  is  contained  completely  within  the  K 
data  samples.  Hence  L  <  K.  A  future  report  will  consider  the  effects  of  signal  segmentation  (i.e.. 
when  the  signal  is  spread  across  a  number  of  sample  blocks  of  length  K).  We  define  an  augmented 
signal  vector  of  length  K  such  that  the  first  L  elements  arc  the  elements  of  s  and  the  remaining 
elements  are  zero.  Let  s'  be  the  resultant  output  vector  after  s  has  been  processed  through  the  GS 
canceller  and  s,j„^.  be  the  resultant  GS  output  augmented  vector.  This  resultant  output  vector  is  then 
input  into  the  matched  filter  of  the  vector  s.  or  equivalently  .  (f  we  '■ct  r'  equal  to  the  response 
of  sj„(,  match  filtered  with  s„„j,  then 


r  '  —  ' 

where  is  defined  as  the  K  x  (2K  -  1)  augmented  ACF  matrix  of 


(12) 


8 


NRL  REPORT  ‘Jl()7 


The  results  and  derivations  to  be  presented  are  the  same  whether  we  use  the  augmented  or  non- 
augmented  notation.  Hence,  we  assume  that  all  vectors  are  augmented  and  drop  the  augmented  desig¬ 
nation.  Note  also  that  we  have  left-justified  .s  in  Actually  .s  could  start  anywhere  in  s,, :  the 

assumption  that  the  s  is  left-justified  is  made  for  convenience's  sake  and  does  not  change  the  results 
of  this  analysis. 

Often  s  is  chosen  so  that  the  matched  filter  response  has  low  sidelobes  ti  e.,  rim)  <  <  mOi  for 
m  y  0).  However,  if  the  desired  signal  is  passed  through  a  GS  canceller  structure,  the  desired  signal 
vector  is  perturbed  and  degradations  occur  in  the  matched  filter  response.  E.xamples  ot  ci'des  that 
have  high  compression  ratios  and  low  sidelobes  are  the  Frank  Code  [12).  Lewis  and  Krctschmei  s 
PI  P4  code  11.1).  and  shift  register  codes  (also  see  Ref.  14).  All  of  these  codes  have  an  ACF  with  all 
sidelobes  well  below  the  matched  response.  Figure  8  shows,  for  example,  the  ACF  of  the  KMI- 
element  Frank  code. 


SAMPLE  NUMBER 


Fig.  8  —  Frank  code  autocorrelation  function  L  =  IIX). 
zero  Doppler  shift  and  no  bandwidth  limitation 


4.  GS  SIMPLIFICATIONS 

This  section  discusses  a  number  of  results  that  significantly  simplify  the  forthcoming  analysis. 
Let  C  be  any  (A  -  1)  x  (A  ~  1)  nonsingular  matrix.  It  is  well  known  |21  that  transforming  the 

auxiliary  input  channels  .v  | . \\  |  by  this  transform  does  not  change  the  transient  or  steadv -state 

noise  power  residue  performance  of  the  SMI  (or  GS).  This  is  because  the  output  from  a  canceller 
with  transformed  auxiliaries  is  identical  to  the  output  from  a  canceller  with  untransformed  auxiliaries 
The  GS  canceller  implementation  is  equivalent  .o  a  matrix  transforming  the  input  channels.  F-or  a  (iS 
canceller,  the  matrix  C  has  the  upper  triangular  matrix  form.  In  the  configuration  illustrated  in  Fig.  4 
the  matrix  transform  C  is  implemented  by  passing  the  input  channels  through  a  GS^.  v  structure  fol¬ 
lowed  by  a  power  equalizer  on  the  output  auxiliary  channels.  The  output  powers  of  the  .M  X  chan 
nels  after  power  equalization  are  equal  to  up,,,,,.  Without  loss  of  generality,  we  can  define  up,,,,  ^  I . 

The  structure  shown  in  Fig.  9  illustrates  that  any  GS  canceller  structure  can  be  divided  into  two 
parts:  a  deterministic  steady  state  front-end  processor,  in  which  the  main  channel  is  decorielated 
from  the  auxiliary  channels,  and  a  stochastic  back-end  processor,  which  is  driven  b\  uncorrelated 
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Fig.  9  —  Residue  equivalent  GS^jv  canceller  using  the  power  equalizer  matrix 


equal  powered  noise  in  each  channel.  The  back-end  processor  is  independent  of  the  input  covariance 
matrix,  and  the  auxiliary  weights  associated  with  the  back  end  processor  go  to  zero  as  AT  —  oo. 
Hence,  the  convergence  properties  of  the  GS  canceller  can  be  studied  by  analyzing  the  convergence 
properties  of  the  back-end  processor.  From  this  point  on,  we  assume  that  the  input  channels  are 
orthogonal  and  of  equal  power. 

A  second  matrix  transform  that  significantly  simplifies  the  forthcoming  analysis  is  now  dis¬ 
cussed.  Let  $  be  any  K  x  K  unitary  matrix,  i.e.  4>'4>  =  !/(,  where  If;  is  the  K  x  K  identity  matrix. 
Let  us  transform  each  input  channel  data  set  x„.  n  =0,  1,2 . N  -  1  by  'i>  such  that 

y„  =  4'x„  .  n  =  0,  I, - N  -  \.  (13) 
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where  y„.  n  =0.  I ,  .  .  .  ,  N  -  1  is  the  resultant  output  data  set.  If  we  input  this  data  set  into  a 
GSx  ,v  canceller  then  [4]  the  estimated  weights  using  the  x„  input  are  identical  to  those  using  the  y„ 
inputs. 

Let 

x„  =  (.r„(l)..v„(2) . x„(K)f  (14) 

be  the  A'-length  input  vector  in  the  nth  channel  and 

z  -  (^(l).c(2) . z(K)V  (15) 


be  the  output  vector  of  a  GSk  i\j.  The  GS kjj  canceller  transforms  the  /f-length  data  vectors  X], 

X2,...,X;v-i  into  an  orthogonal  set  of  AT-length  data  vectors  Z|,  Zt . z.v-i-  In  fact,  it  was 

shown  in  Ref.  15  that 


z  =  Gxo, 

where  G  is  the  GS  complementary  projection  matrix  given  by 


Ziz'i  Z-)Z'-> 

G  =  Ik-^ 

z'lZ]  Z2Z2 


Z/V  -  I  zfv  -  I 

Zv-lZv-l 


(16) 


(17) 


It  was  also  shown  that  G  can  be  written  as 


G  =  4>'  Ayv-i***.  (18) 

where  4>  is  a  A  x  A  unitary  matrix  that  is  a  function  of  x^,  n  =  1,  2,  .  . .  ,  A-1  and  A,v-i  is  a 
K  X  K  diagonal  matrix  with  the  first  N  -  1  diagonal  elements  equal  to  zero  and  the  other  diagonal 
elements  equal  to  one.  Thus 


z  =  4>'Aa,_,  ^>Xo. 

We  can  write  the  A-length  input  vector  in  the  main  channel  as 


(19) 


Xm  =  s  +  Xo, 


(20) 


where  s  denotes  the  A-length  vector  desired  signal  subcode  and  Xq  the  A-length  noise  vector.  Because 
of  linearity,  the  GSj^-jy  canceller  can  be  decomposed  as  shown  in  Fig.  10.  Here  the  left  hand  GS;<  v 
canceller  has  only  the  desired  signal  in  the  main  channel  and  the  right-hand  GS^-  ;v  has  only  Xo  in  the 
main  channel.  Note  that  the  interior  weights  of  the  GS^jv'  are  not  identical  because  of  the  different 
main  channel  in  each  (actually  only  the  weights  in  the  main  channel  differ).  Hence,  the  left-hand 
GSk  1^  output  contains  the  perturbed  desired  signal  output  and  the  right-hand  GS;^-  v  output  contains 
the  output  noise  residue.  For  the  forthcoming  analysis,  we  consider  only  the  left  hand  output  that 
contains  the  desired  signal. 
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Fig.  10  —  Decomposition  of  signal  and  noise  for  concurrent  processing 


Figure  9  shows  that  the  desired  signal  passes  unperturbed  through  the  front-end  processor 
GSoo  /v-  This  is  because  the  desired  signal  vector  is  a  constant  and.  hence,  is  uncorrelated  with  the 
auxiliary  channels.  Thus  the  main  channel  weights  in  the  GSoo,/v  canceller  are  zero.  As  a  result,  we 
need  only  consider  the  effects  on  the  desired  signal  by  the  back-end  processor  where  all  the  auxiliary 
channels  are  equipowered  and  independent  random  variables. 

Finally,  it  was  shown  in  Refs.  3  and  4  that  if  Sq  and  s,v-i  are  the  input  and  output  signal  vec¬ 
tors  respectively,  of  a  GSy^  and  the  auxiliary  inputs  satisfy  assumptions  1  through  5  (Section  2). 
then  T]  =  Sv-iS^_i /SqSq  has  the  following  probability  density  function  (PDF); 


Piv)  = 


(K-\y. 


(K-NV.  (A -2)! 


0  <  77<  1. 


(21) 


From  this  it  is  straightforward  to  show  that  for  K>N, 


fjS^v-lS^-l  1 


i  - 


A-1 

K 


SoSo 


(22) 


and 


^1  I  S/V-|S^r_l  I  ‘1  - 


1 


2(A-I)  ^  A(A-l) 
K  K{K  +  \) 


SoS(i  I  ' 


(23) 


and  El  I  denotes  expectation. 
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5.  SINGLE  AUXILIARY  CROSS-CORRELATIONS 

This  section  presents  derivations  for  two  cross-correlations  of  outputs  of  a  2-input  GS  canceller. 
These  expressions  are  fundamental  in  solving  for  the  average  pulse  compression  sidelobe  increase 
when  N  >  2. 

Consider  two  separate  augmented  signal  vectors  Uq  and  Vq.  We  assume  that  the  desired  signal 
vectors,  each  of  length  L,  are  completely  contained  in  Uq  and  Vq,  respectively.  The  augmented  signal 
vectors  have  length  K.  Thus  L  <  K.  We  input  Uq  and  Vq  into  separate  GSfc  i  cancellers  with  a  com¬ 
mon  auxiliary  vector  x  (Fig.  II).  The  K  samples  (or  elements)  of  x  are  random  variables  that  satisfy 
assumptions  1  through  5  (Section  2). 

The  output  signal  vectors  of  the  respective  GS^  2  (P'g-  U)  denoted  by  Uj  and  V|.  We  wish 
to  find  £1  |u'i  V]  1^}  and  ElUui  ||^|lvi  ||^j  where  for  any  vector  c,  ||c|p  =  c'c. 


Fig.  1 1  —  Parallel  2-input  GS  cancellers 


From  the  previous  section,  we  know  that  any  KxK  unitary  matrix  transformation  of  each  input 
vector  into  a  GS  canceller  does  not  change  the  GS  weights  that  are  calculated  thereafter.  Further¬ 
more.  it  is  elementary  to  show  that  if  U|  and  v,  are  the  respective  outputs  of  the  unitary  matrix 
transformed  canceller,  then 

E\  |u',v,  1^1  =  Ellu'iv,  1^1,  (24) 

and 

f(l|uif  l|vi||^l  =  £l||ii|||^  ||v,|p).  (25) 

We  use  this  relationship  to  simplify  the  input  signal  vectors  Uq  and  Vq.  It  can  be  shown  [16) 
that  3i  K  X  K  unitary  matrix  ^  exists,  such  that 

Uo  =  <I>  Uq  =  lluoild,  0,  0 . 0)^  (26) 

and 

Vo  =  4»Vo  =  IIvoIKp,  IpI^  ,  0,  0,  ,  0)^,  (27) 
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where 


^  lluoll  llvjl  ’ 

and  we  have  defined  Uo  and  Vq  as  the  unitary  matrix  transformed  input  vectors.  Set 

y  =  ^»x.  (29) 

Note  that  because  x  is  a  normalized  ^-length  multivariate  complex  circular  Ga..i3i.*r  vector, 
y  =  O’l,  >’2.  •  •  ■  '  y'Kf  is  also. 

The  ^-length  output  vectors  of  each  GSf^  2  are  given  by 


~  '  y  '*0  ,  yy*  ~ 

Ui  =  Uo - y  =  Ik  -  ^  uo 

y  y  y  y  ) 


y  Vq  yy' 

Vi  =  Vo - —  y  =  Ik  -  ^  Vo, 

y  y  y  y  J 


where  Ik  is  the  K  x  K  identity  matrix.  We  can  show  that 


U|V|  =  lluoll  llvoll  p  1  -  ^ 

I  yy 


1  -  Ip 


ju.p  =  Kr  1  -  ^ 

y  y 


'v,f  =  llvoll'  I  - 


2  „  |2 


p*y,  +  Vl  _  I  p  I  ^  y2 


It  is  shown  in  the  Appendix  that 


,2.  ,  /  ,2  f  .  2  ,  1  ,  liuoll^  llvoll' 

£||„,v,|-|=£||u,v,!  1=  lu„v„|  +  ,) 


IuqVoI^  .  „  .  2  .  1 


E\\\u,f  tiv.ri  =  £iiii,ir^  iiv.in  =  +  iKip  iivor  i  -  ^  + 


K{K  +  1) 


K  K(K+\) 
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6.  MULTIPLE  AUXILIARY  CROSS-CORRELATIONS 


In  this  section,  we  derive  expressions  for  two  cross-correlations  of  outputs  of  an  /V-input  GS 
canceller.  The  input  signal  vectors  Uq  and  Vq  are  as  defined  in  the  previous  section.  We  input  Uq 
and  Vq  to  separate  GS^^-jv  cancellers  with  common  auxiliary  vectors  Xi,  Xt,  ...  ,  x,v_i  as  illustrated 

in  Fig.  12.  The  K  samples  (or  elements)  of  each  auxiliary  vector  X| . x^-i’  ^^e  r.v.s  that 

satisfy  assumptions  1  through  5  (Section  2).  Let  the  signal  output  at  each  level  of  the  respective 
GS^jv  be  denoted  by  u„  and  v„,  «  =  0,  1,  . . .  ,  N  —  1.  We  wish  to  find  £[  |  u^_iVv-i  |  "I  and 
|p||v^_i  ll^j.  It  will  be  found  that  this  can  be  done  recursively. 

Uq  Xi  X2  Xm_1  Vq  Xg  ^N-1 


To  this  end,  as  in  the  previous  section,  we  multiply  Uq  and  Vq  by  a  unitary'  matrix  transforma¬ 
tion  <I>,  such  that  the  new  input  vectors  Uo  and  Vq  are  given  by  Eqs.  (26)  and  (27).  It  is  straightfor¬ 
ward  to  show  that 


E\ 


E\ 


l|2i  _ 


=  f|||u„ 


r.  II  2  II  ~ 


(37) 


for  n  =  0,  1 ,  ....  A  -  1 .  In  addition,  each  of  the  auxiliary  input  data  vectors  is  multiplied  by  <I> 
such  that  x„  =  4>x„,  n  =  1 .  2,  . .  .  ,  -  1 .  Note  that  each  of  the  transformed  auxiliary  vectors  also 

satisfies  assumptions  1  through  5  and  is  statistically  identical  to  x„,  n  =  1 ,  2,  .  . .  ,  N  -  1 . 

We  redraw  the  configurations  seen  in  Fig.  12  into  the  equivalent  configurations  seen  in  Fig.  13. 
Here  we  show  for  each  of  the  original  cancellers  decomposed  into  two  parallel  can¬ 

cellers  followed  by  a  single  GS^  2  canceller.  The  /C-length  noise  vector  into  the  GS^:,  2  cancellers  is 
denoted  by  y. 

As  stated  in  Section  4,  we  can  write 


U;V-2  —  ^0  ^N-2  ^0  Uq* 

(38) 

^N-2  =  ^0  -2  ‘J’O  Tq. 

(39) 

y  =  4»o  A^_2  ^0  *1' 

(40) 
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Fig.  13  —  Equivalent  representation  of  Fig.  12 

where  a  K  x  K  unitary  matrix  that  is  a  function  of  X2,  x  ,  . . .  ,  Xn-]  and  A;v-2  a  /ii  x 
diagonal  matrix  with  the  first  N  -  2  diagonal  elements  equal  to  zero.  All  other  diagonal  elements  are 
equal  to  one. 

If  we  premultiply  U/v_2,  V;v-2>  y  the  unitary  matrix  the  resultants  of  the  cross¬ 
correlations  we  are  seeking  do  not  change.  However,  the  number  of  independent  and  identically  dis¬ 
tributed  (i.i.d.)  elements  in  the  auxiliary  vector  is  reduced  by  A  -  2  as  indicated  by  Eq.  (40).  The 
first  N  -  2  elements  of  4>y,  4>U/v-i,  and  <I>v^_2  are  zero.  Hence,  it  follows  using  the  methodology 
given  in  the  previous  section  for  finding  cross-correlations  for  an  arbitrary  number  of  i.i.d.  samples  in 
the  auxiliary  channel  of  a  2-input  GS  that  the  expectations  of  these  cross-correlations  conditioned  on 
X2,  ■  ■  ■  ,  -Xv-i  are  given  by 


£^1  I  Uv  -  I  V/V  - 1  I  ^  I  X2 ,  .  .  . ,  X;v  _  1  ! 


(41) 


UAr-2''/V-2 


1  - 


-t- 


1 


K  -N  +  2  (K  -N  +  2)(K  -  A  +  3) 


-I- 


llUA/-2lrl|Vv 


{K-N  +  2){K  -  A  +  3) 


and 


£illUAf-lll  1*2'  •  •  •'  X/V-l: 


(42) 


UV-2V/V-2  1  ^ 


(A  -  A  +  2){K  -  A  +  2) 


+  i|U/V-2ll^l|Vv- 


+ 


1 


(K  -N  +  2)  (K  N  +  2)(K  -  N  +  3) 
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Wc  rewrite  Eqs.  (48-50)  in  vector  form  by  defining  e„  =  (a„,  as 


A 

1  UqVo  I  ^ 

Cn  +1  ~ 

6„  ,  I.C.;  €o  = 

K  IP  11  Vo  f 

Thus  solving  for  f/vr_i : 


N-2 

b„' 

- 

'  !u{,Vo|2  ■ 

CA-1  - 

n 

n  =0 

b„  a„ 

KlPlIvoll^ 

(51) 


(52) 


The  desired  cross  correlations  are  given  by  the  elements  of  eyv-i  where  £j  |  u5v-iv,v_i  |‘t  and 
Ejllu/v-i  ||“||v/v-i  11^1  are  given  by  the  first  and  second  elements,  respectively. 

Detine 

T/i ,,(/:,  /V)  a,2(K,  n) 

^  =  A2i(K,  N)  A22iK,  N) 

Then 

E\\u'n-\Vn-\  1^1  =  ^ii(/^,  N)  |u(,Vo|^  +  Ai2(K,  N)  l|iiolh|vof  (54) 

and 

£'ll|uA-if  l|vv-i|pl  =  A2i(K,  N)  |u{)Vo|^  -t-  A22iK,  N)  l|uof  ||vof.  (55) 

We  note  that  because  of  the  symmetries  in  form  of  the  2  x  2  matrix  A  defined  in  Eq.  (53).  we 
can  show  that  A2i(K,  N)  =  Ai2(K,  N).  In  fact,  if  we  evaluate  Eqs.  (54)  and  (55)  at  Uq  =  vq,  we 
find  that 


N-2 

bn 

=  n 

n  =0 

bn 

On 

(53) 


Au(K.  N)  =  A22{K,  N) 


(56) 


and  using  Eq.  (23) 


Ai^iK,  N)  +  A^2(K,  N)  =  1 


2iN  -  1)  N(N  -  1) 

K  K{K  +  \)' 


(57) 


7.  SIDELOBE  DEGRADATION  DERIVATION 

We  now  derive  an  exact  expression  for  the  average  adaptive  pulse  compressed  sidelobe  by  using 
the  result  of  the  previous  section.  The  desired  input  signal  vector  is  assumed  not  to  be  segmented 
(i.e.,  it  is  contained  entirely  in  the  augmented  signal  vector  denoted  by  Sq).  Each  auxiliary  channel 
has  K  samples  and  these  samples  satisfy  assumptions  1  through  5.  Also,  the  desired  signal  vector  is 
normalized  so  that  ||so||^  =  1 
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Figure  14  shows  the  A^-input  canceller  followed  by  the  matched  filter.  Let  s^.  be  defined  as  any 
column  of  the  augmented  ACF  matrix.  Then,  an  expression  representing  an  output  r  of  the  matched 
filter  (match  point  or  sidelobes)  after  cancellation  is  given  by 

r  —  s'.S;v_i  (58) 

where  S/^-^  is  the  ^-length  output  signal  vector  of  the  GSkjj  canceller. 

Sq  X.,  X2  ^N-1 


We  define  the  average  adaptive  pulse  compression  sidelobe  level  associated  with  s,,  as 
which  is  given  by 


£1  |s^s.v_i  |-i 

£|  I  S/v-iS/v_i  I  “I 


(59) 


As  A"  —  00.  note  that  s^i-i  So  and  5£„ goes  to  the  quiescent  sidelobe  level  of  s, .  This  equals 
Isl  sol"  which  is  given  as  S£u We  normalized  by  the  adaptive  match  point  response 
£1  1  S/v-iS/v-t  |“L  just  as  ]  sJ.Sq  r  would  be  normalized  with  the  quiescent  match  point  response 

I  SqSo  I  “  . 

The  denominator  of  Eq.  (59)  was  evaluated  in  Refs.  3  and  4  and  in  Section  4  and  is  given  by 


£j|sV-,SA,_,  =  1  - 


2(N  -  1) 
K 


N{N  -  1) 
K(K  +  1) 


(K  -  N  +  l)(K  -  N  +  2} 
K(K  +  1) 


(60) 
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The  numerator  of  Eq.  (59)  is  evaluated  by  using  the  results  of  the  previous  sections.  As  stated 
in  Section  4,  if  Sq  and  s,.  are  the  main  channel  inputs  to  two  distinct  GSx.n  cancellers  with  identical 
auxiliary  inputs,  then  the  outputs  of  these  two  cancellers  can  be  written  as 


S^_l  —  4>'Ayv-l4’  Sq, 

^c.N-l  =  So 

where  is  a  A  x  K  unitary  matrix.  Using  Eqs.  (61)  and  (62),  it  can  be  shown  that 

Sr.V-lSN-l  =  S^S/v-l. 

Thus 

E{ \s',SN-y  =  £■{  |s^./v-iSv-i  1^1. 


(61) 

(62) 


(63) 


(64) 


and  the  analysis  of  Section  6  can  be  applied  with  Uq  =  s<.  and  Vq  =  Sq. 

Using  Eq.  (54)  and  the  fact  that  ||so||“  =  1,  it  follows  that 

ftls's^-i  1^!  =  A^^{K,  N)  |s;.So|2  +  A,2{K,  N)  11s,11'. 


Substituting  Eq.  (60)  and  (65)  into  Eq.  (59)  results  in 
K{K  +  l)/ln(A:,  N) 


SL. 


iK  -  N  +  \){K  -  N  +  2) 


SL„ + 


\(K  +  \)A^2(K.  N) 

{K  -  N  +  \)(K  -  N  +  2) 


(65) 


(66) 


We  define 


<2(A,  N)  = 


K{K  +  \)Au{K,  N) 

{K  -  N  +  \)(K  -  N  +  2) 


to  be  the  quiescent  sidelobe  factor  and 


ASLJK,  N) 


K(K  +  \)  A  N) 

(A  -  A  +  1)(A:  -  N  +  2) 


(67) 


(68) 


to  be  the  adaptive  sidelobe  perturbation.  Note  ASL„{K,  N)  is  constant  over  all  sidelobes  of  the 
compressed  pulse. 

8.  RESULTS 

For  the  results  to  be  presented,  it  is  convenient  to  set  =  N  -  \  where  is  the  number 
of  auxiliary  channels.  In  addition,  the  number  of  independent  samples  per  channel  K  is  set  equal  to 
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an  integer  multiple  of  Let  K  =  MN^,u.x  Qo(M,  A^„u.,)  =  Q(K,  N)  where  Q{K.  N)  is  defined 

by  Eq.  (68).  To  achieve  good  signal  detection  and  noise  cancellation,  for  most  practical  cases 
A/  >  2,  [2-5,17].  Figure  15  presents  the  quiescent  sidelobe  factor  Qo(M.  plotted  in  dB  vs 
and  M.  Note  for  all  A/  >  2  and  that  is  less  than  and  approximately  equal  to  one. 

Hence  the  approximation  that  Q^)(M.  /V)  =  1  is  valid  in  mo.st  cases.  Reexamining  Eq.  (67)  indicates 
that  for  M  >  1.  the  adaptive  sidelobe  level  is  equal  to  the  quiescent  sidelobe  level  plus  the  adaptoe 
sidelobe  perturbation  term  given  by  Eq.  (69)  which  wi  denote  as  ^.SLjK.S)  ||s, 


NO.  OF  AUXILIARIES 


Tie.  I-S  UuieMoni  MJelubc  huli't  .V .  .iiui  \ 

Assume  max  !|s,  H'  =  1.  Figure  16  plots  max  A.V/.„  vs  and  A/.  These  curses  establish  a 
performance  floor  for  the  peak  adaptive  pulse  compression  sidelobes.  F(>r  a  given  and  K  (or 
A/A',„^, ).  the  adaptive  compressed  sidelobe  level  on  the  average  is  greater  than  the  A.5/.„  calculated  lor 
and  K.  Th*'  figure  shows  that  max  A5Z,„  decreases  monotonically  with  M  and  .V.,,,.  Hence, 
two  ways  of  decreasing  the  performance  floor  are  to  increase  either  the  number  of  independent  input 
samples  or  the  number  of  auxiliars'  channels.  We  also  observe  that  more  than  5  indcqx'ndent 

samples  are  necessary  to  achieve  a  modest  adaptive  compres.scd  sidelobe  level  of  -25  dB  for  leuer 
than  20  auxiliaries. 

It  is  significant  to  note  that  even  had  the  waveform  been  pulse-compressed  intii  a  single  range 
bin  before  the  cancellation  process,  the  maximum  range  sidelobe  level  after  cancellation  w<Hjld  be  the 
same  as  if  the  waveform  had  been  compressed  after  the  canceller  To  see  this,  assume 

s  =  (1.0.0 . 0)^.  This  input  vector  corresponds  to  the  output  of  a  perfect  (uilse  compressor,  i  c  . 

all  of  the  signal  energy  is  in  cne  range  ceil  (we  are  assuming  a  sampling  rate  equal  to  the  uaseform 
bandwidth).  If  this  waveform  is  input  to  the  cancellei.  processed,  and  matched  filtered  with 

s  =  (1.0 . 0)^.  then  Eq.  (67)  indicates  that  this  input  signal  vector  has  maximum  range  suiclobc 

levels  equal  to  ISLjK.N).  Hence,  the  same  perturbed  sidelobe  level  results  regardless  ot  whether 
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auxiliaries 


Fig.  16  -  Maximum  adaptive  sidelobe  penurbations  vs  and  M. 
signal  completely  contained  in  K  input  samples 


pulse  compression  (Kcurs  before  or  after  cancellation.  Again,  note  that  we  are  assuming  the 
uncompressed  pulse  is  com.pletely  contained  in  the  K  samples  used  to  calculate  the  canceller  weights. 
Thus,  it  would  seem  desirable  if  L  «  K  to  place  the  pulse  compression  after  the  canceller  because 
of  the  previously  cited  disadvantages  (see  Section  1)  of  having  it  before  the  canceller.  If  the 
uncompressed  waveform  extends  over  a  number  of  batches  of  K  samples,  it  might  be  more  desirable 
to  do  pulse  compression  first  to  prevent  the  compressed  pulse's  canceller-induced  higher  sidelobes 
from  extending  over  these  batches  A  future  study  will  address  this  issue. 

The  above  analysis  of  pulse  compression  and  canceller  interactions  can  also  be  applied  to  quanti¬ 
fying  the  canceller  degradation  caused  by  the  presence  of  a  desired  signal  in  the  samples  used  to  cal¬ 
culate  the  adaptive  canceller  weights.  If  the  desired  signal  has  the  power  o;  after  pulse  compression, 
then  the  average  power  residue  from  the  signal  in  the  AT  -  1  range  bins  not  containing  the  signal  can 
be  shown  from  our  analysis  to  equal  at  most  a;  ASL„(K.N)  plus  possibly  the  signal  power  resulting 
from  the  quiescent  compressed  sidelobes.  Now  a;  ASL„{K.N)  might  be  much  greater  than  the  quies¬ 
cent  output  noise  power  level  of  the  canceller.  However.  ASL„{K.N)  is  a  function  of  N  and  K. 
and  with  increasing  K  can  be  made  as  small  as  desired.  Hence,  by  using  the  values  of  a;.  N.  and 
a;„in  in  conjunction  with  Fig.  16.  one  can  fitid  the  iiuiiidci  of  independent  samples  K  needed  such  that 
a;  ASL^iK.N)  is  within  3  dB  of  a^,n. 

We  note  that  in  some  canceller  schemes  the  adaptive  weights  are  computed  from  the  K/2  sam¬ 
ples  before  and  after  a  chosen  sample.  Hence  the  chosen  sample  itself  is  not  used  in  the  weight 
evaluation.  This  technique's  advantage  is  that  if  the  desired  signal  is  kKalized  to  a  single  range  sam¬ 
ple  the  undesirable  effect  of  cancelling  the  desired  signal  does  not  (Kcur  13.4).  However,  a  noise 
power  residue  from  the  desired  signal  is  present  in  the  canceller  output  samples  preceding  and  follow¬ 
ing  the  chosen  sample  because  the  signal  is  present  in  the  weight  calculation  for  the‘e  samples.  This 
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level  occurs  in  the  K H  samples  before  and  after  the  chosen  sample.  The  ma.Kimum  level  of  the 
signal-induced  canceller  residue  is  equal  to  o"  ASL,,(K,N).  Again.  K  can  be  chosen  large  enough  to 
be  within  3  dB  of 

One  final  note:  for  some  applications,  the  matched  filter  is  replaced  by  a  filtering  scheme 
whereby  the  range  sidelobes  are  reduced  at  the  e.xpense  of  signal  gain  at  the  match  point.  However, 
the  results  derived  in  this  report  are  also  valid  if  any  filter  other  than  the  matched  filter  s,,  is  used. 
We  could  replace  the  So  seen  in  the  "matched  filter"  block  in  Fig.  14  with  a  general  weighting  func¬ 
tion  given  by  the  T-length  vector,  a  with  elements  Oo.  0| . lo  oor  analysis,  we  would 

replace  the  S  matrix  defined  by  Eq.  (10)  with  an  A  matrix  whose  elements  are  given  by  replacing  the 
,v  s  with  a  s  in  Eq.  (10).  The  vector  s,  then  would  be  taken  to  be  any  column  in  A  and  the  anahsis 
follows  as  given. 

9.  SLiM.MARV 

Performance  results  for  the  sidelobe  level  of  a  compressed  pulse  that  has  been  preproces^cu 
tnrough  an  adaptive  canceller  have  been  obtained.  The  adaptive  canceller  is  implemented  b\  u^ing 
the  Sampled  Matrix  Inversion  (SMI)  algorithm.  Because  of  finite  sampling,  the  quiescent  comprc'"'ed 
pulse  sidelobe  levels  are  degraded  because  of  the  preprocessing  of  the  main  channel  input  data  '•tream 
(i.e.,  the  uncompressed  pulse)  through  an  adaptive  canceller.  An  exact  expression  was  derived  tor 
the  average  adaptive  pulse  compression  sidelobe  level  that  was  hvund  to  be  a  function  of  the  nutiiber 
of  auxiliary  input  canceller  channels  the  number  of  independent  samples  per  channel  K.  and  the 
given  input  code  of  length  L.  It  was  shown  that  K  /  can  be  significantly  greater  than  one  to 
retain  sidelobes  that  are  close  to  the  original  quiescent  sidelobe  level  (with  no  adaptive  canceller) 
Also,  it  was  shown  that  the  maximum  sidelobe  level  degradation  is  independent  of  whether  puNe 
compression  occurs  before  or  after  the  adaptive  canceller  if  the  unctmipressed  pulse  is  completelv 
contained  within  the  K  samples  that  are  used  to  calculate  (he  canceller  weights.  It  was  further  shown 
that  this  same  analysis  can  be  used  to  predict  the  noise  power  level  at  the  canceller  output,  that  is 
induced  by  having  the  desired  signal  present  in  the  canceller  weight  calculation. 
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Appendix 

DERIVATION  OF  EQS.  (35)-(36) 


Starting  with  Eq.  (32),  we  first  derive  Eq.  (35).  Now 
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Note  that  the  expected  values  of  the  cross-terms  resulting  from  magnitude-squaring  the  right  side  of 
Eq.  (32)  are  zero  because  .V2  is  zero  mean.  Now  for  p  =  1,  we  can  show  using  Eqs.  (Al)  and  (23) 
that 
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From  Eq.  (A2),  it  can  be  shown  that 
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Because  y|,  ..,  }’/(  are  identically  distributed  r.v.s.  the  terms  in  the  first  summation  are  all 

equal  and  the  terms  in  the  double  summation  are  equal.  Let 
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Substituting  Eqs.  (A3)  and  (A5)  into  Eq.  (A4)  results  in 
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Thus, 
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Equation  (35)  follows  directly  from  Eq.  (A8). 

We  now  derive  Eq.  (36).  From  Eqs.  (33)  and  (34) 
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Now  the  expectations  of  the  individual  terms  of  Eq.  (A9)  can  be  evaluated  by  using  some  of  the  pre¬ 
vious  results. 
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Thus, 
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Equation  (36)  follows  from  Eq.  (A  13). 
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